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ABSTRACT 

We present a method for the identification of heavily absorbed AGN (iVn> 10 23 cm~ 2 ) 
from X-ray photometric data. We do this using a set of XMM-Newton reference spectra 
of local galaxies for which we have accurate A"h information, as described in Brightman 
&: Nandra. The technique uses two rest-frame hardness ratios which are optimised 
for this task, which we designate HR1 (2-4/1-2 keV) and HR2 (4-16/2-4 keV). The 
selection method exploits the fact that while obscured AGN appear hard in HR2 due 
to absorption of the instrinsic source flux below ~ 4 keV, they appear soft in HR1 due 
to excess emission originating from scattered source light, thermal emission, or host 
galaxy emission. Such emission is ubiquitous in low redshift samples. The technique 
offers a very simple and straightfoward way of estimating the fraction of obscured 
AGN in samples with relatively low signal-to-noise ratio in the X-ray band. We apply 
this technique to a moderate redshift (z ~ 1) sample of AGN from the Chandra Deep 
Field North, finding that 61% of this sample has Ah> 10 23 cm -2 . A clear 
and robust conclusion from our analysis, is that in deep surveys the vast majority of 
sources do not show hardness ratios consistent with a simple absorbed power law. The 
ubiquity of complex spectra in turn shows that simple hardness ratio analysis will not 
yield reliable obscuration estimates, justifying the more complex colour-colour analysis 
described in this paper. While this method does very well at separating sources 
with Ah> 10 23 cm 2 from sources with lower Ah, only X-ray spectroscopy 
can identify Compton thick sources, through the detection of the Fe Kq 
line. This will be made possible with the high throughput X-ray spectral 
capabilities of ATHENA 
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1 INTRODUCTION 

The nature and amount of obscuration in active galactic nu- 
clei (AGN) is an important issue for a range of astrophysical 
questions, such as unification of AGN (eg. |Antonucci||1993[) 
and the role of AGN in galaxy evolution ( |Hopkins et aL] 
20061. Heavily absorbed AGN can be identified using their 



2-10 k eV band is significant once th e Ah reaches a few 10 



X-ray spectra, via direct measurement of the line of sight 
column density. At the highest column densities (Ah > 10 24 
cm" 2 , where even the direct X-rays are suppressed, obscured 
AGN can still be identified via the presence of X-ray reflec- 
tion, with the spectrum consisting of a flat X-ray continuum 



accompanied by a high equivalent width Fe Kcv line (Matt, 
Brandt fe Fabian|1996 \ . The measurement of column densi- 



ties or identification of reflection features requires good qual- 
ity X-ray spectra with enough counts to carry out spectral 
fitting and to reject alternative hypotheses. This is challeng- 
ing, even in the local universe, as the flux suppression in the 



E-mail: mbright@mpe.mpg.de 



cm (Brightman & Nandra 2011al. Many alternative meth- 



ods of identifying heavily absorbed sources, where the X-ray 
spectrum itself is not good enough to do so, have been pre- 
sented in the literature. These often use the ratio of the ob- 
served X-ray flux to the flux at some wavelength at which the 
effect of the obscuration is assumed to be small, and which 
should therefore be representative of the AGN bolometric 
flux. Examples of these assumed "isotropic indicators" are 



the MIR, [Om] line or [Nev] line ( Diamond-Stanic et al 



2009}|LaMassa et al.|20l"o"||Gilli et al.|2010| >. Obscured AGN 
may be found via weak observed X-ray emission with respect 
to the isotropic indicator, with the X-ray flux suppression 
attributed to absorption. A popular technique which uses 



this method was presented by Bassani et al. ( 1999 1 and uses 



Fx/Fiqm] in conjunction with the Fe Kq line equivalent 
width as proxies for A^h, though -Fx/-F[oni] has henceforth 
often been used alone. The use of isotropic indicators are 
complicated when the isotropic indicator has contributions 
from non nuclear sources (e.g. star-formation processes con- 
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tribute to both MIR and [O ill] fluxes), or where the source is 
intrinsically X-ray weak. Furthermore, techniques have even 
been developed that use data entirely outside of the X-ray 



band to detect X-ray obscuration (e.g. Daddi et al. 2007 



Fiore et al.|2008 1. Arguably though, the best place to judge 



X-ray absorption, is in the X-ray band alone. This is some- 
times done using 'hardness ratios', which is a model indepen- 
dent technique, and utilises X-ray photometry. Traditionally, 
hardness ratios are calculated using observed counts, often 
in the standard soft (0.5-2.0 keV) and hard (2.0-10.0 keV) 
bands. Sources with large hardness ratio are assumed ab- 
sorbed, while sources with low ratios assumed unabsorbed. 
However this method is very crude and may not exclusively 
separate the obscured sources from the unobscured sources if 
spectral complexity is present. At other wavelengths, colour- 
colour selection is an effective method of identifying sources 
populations based on photometry. Here, we seek a novel ap- 
proach to the identification of X-ray absorption from X-ray 
photometric data, using colour-colour selection, formed from 
X-ray hardness ratios. While work using a combination 
of hardness ratios to infer absorption has been pre- 
sented previously, for example results by Delia Ceca 



et al. (1999) using ASCA data, we seek to improve on 



these using knowledge from local samples with data of high 
spectral quality. This better understanding of hardness ra- 
tios can then be applied in cases with much poorer data qual- 
ity, such as in deep surveys of the distant Universe, and help 
elucidate the observed relationships seen between obscura- 
tion, luminosity and redshift (eg. |Ueda et al.|2003 Hasinger 



2008; Brig htman fe Nandra||2011b[ Brightman fe Ueda, in 
preparation 



We use a physically and empirically motivated choice of 
X-ray bands for our colour-colour plots, with consideration 
for the Fe Ka line and the expected energy of the absorption 
turn over in heavily obscured AGN. In Fig [l] we present a 
selection of our reference XMM-Newton model spectra, over 
a range in Ah values. It can be seen from these that for 
iV H > 10 23 cm" 2 there is little direct AGN flux observed be- 
low 2 keV, and that emission below 2 keV is dominated by a 
soft-excess component, which can be attributed to scattered 
AGN emission, or emission from the host galaxy. As the line 
of sight Ah approaches 10 24 cm" 2 , this galactic/scattered 
emission can dominate even up to 4 keV in the rest frame, 
and the spectral characteristics of the heavily obscured AGN 
is seen only above 4 keV. This can be a reflection dominated 
spectrum, characterised by its flat shape, or a transmitted 
spectrum with high measured TVh- Both of these are usually 
accompanied with iron lines of high EW. Based on these 
physical motivations, we choose our energy band selection 
to be 1<E<2 keV (1), 2<E<4 keV (2) and 4<E<16 keV 
(3). To analyse the colours, we then extract model photon 
fluxes from our set of reference spectra. While the 4-16 
keV data include an extrapolation out of the XMM- 
Newton band, from 10-16 keV, the spectra are well 
modelled and the iVn well denned due to high qual- 
ity data below 10 keV, making uncertainties in the 
extrapolation minimal. We use fluxes rather than count 
data so that the selection technique is independent of the 
telescope or instrument used. We calculate hardness ratios 
based on these fluxes, where HR1= Z^ and ,V\^ ^t an j]l and 

' F(band2) + F{ banal) 

„n _ F(band3)-F(band2) 



2 X-RAY DATA 

For the task of finding an effective colour-colour selection 
technique for obscured AGN, we utilise a set of 126 XMM- 
Newton spectra of local galaxies, for which we have carried 
out a detailed spectral analysis of in |Brightman fe Nandra] 
( |2011a| BN11). The parent sample of this set was the ex- 



tended IRAS 12 micron galaxy sample of Rush et al. ( 1993 1 



The mid-IR selection of this sample should ensure relatively 
little bias against obscured objects ( |Spinoglio fe Malkan| 
19891. The X-ray spectra of these galaxies were fitted in 



a systematic way, starting with a simple power-law, and 
adding additional model components, such as absorption or 
reflection when statistically required. BN11 also applied new 
Monte-Carlo based models for modelling heavily absorbed 
sources. Compton thick sources were identified either with 
a direct measurement of the line of sight A^h, or from the 
presence of a reflection dominated spectrum, often accom- 
panied by strong Fe Ka emission. As such, we have accurate 
]Vh information for these galaxies from XMM-Newton. As 
the sample is local (z< 0.1), we treat the X-ray spectral 
data as rest-frame. The sample consists of all galaxy types, 
from Seyfert galaxies to Hn and LINER galaxies, with a 
wide range of Lx, A^h and a high fraction of absorbed and 
Compton thick sources. For a more detailed description of 
the X-ray spectral analysis, the X-ray spectral properties of 
the sample, and the new Monte-Carlo models, see BN11. 
For the optical activity classification and 12 urn data, see 



Brightman & Nandra (2011b I 



3 COLOUR-COLOUR SELECTION FOR 
HEAVILY OBSCURED AGN 

In Fig.[2]we present a colour-colour plot which uses our hard- 
ness ratios using our reference spectra data. On it, we also 
show the Ah of each source, represented by colour. There 
is a cluster of sources at HRl~-0.25 and HR2~0.0. These 
are objects which show little absorption above 1 keV. Both 
ratios harden as the A^h increases until ~ 10 23 cm" 2 . At this 
point, the soft-excess emission begins to dominate below 4 
keV and hence HR1 softens, while HR2 remains hard. It can 
be seen that the heavily absorbed sources (A H > 10 23 cm" 2 ) 
are separated well from the others using these hardness ra- 
tios in combination, being hard in HR2, while soft in HR1. 
We therefore plot two proposed selection criteria. First (1), 
the dotted line selects almost all (38/41=93%) heavily ab- 
sorbed sources, with only a few (3/41=7%) contaminants 
with a lower Ah- Secondly (2), the solid wedge selects all 
(16/16=100%) Compton thick sources. This region is con- 
taminated by 20/36 (56%) sources with A H < 10 24 cm" 2 , 
though only by 3/36 (8%) sources with A H < 10 23 cm" 2 . 
The selection criteria are as follows. For Ah> 10 23 cm" 2 , 

(1) HR2 > 0.62 x HR1 + 0.38 (1) 
and for Compton thick AGN 

(2) HR2 > 0.25 AND HR2 < 2.28 x HR1 + 0.36 (2) 

Calculation of the rest-frame hardness ratios is in prac- 
tice very straightforward. For example, counts or count 
rates can be converted to flux to be used in this selection 
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Figure 2. Our colour-colour plot using the hardness ratios HR1 (2-4/1-2 keV) and HR2 (4-16/2-4 keV), calculated from model photon 
fluxes extracted from our set of 126 XMM-Newton reference spectra. The colour of each data point corresponds to the value of the 
measured jVjj in that spectrum. The dotted line and the solid wedge are proposed selection regions for heavily obscured AGN, as 
described in the text. 



technique using the HEASARC Portable Interactive Mulit- 
Mission Software (pimmsQ This requires the assumption of 
a model. For deep X-ray fields, the use of a F = 1.4 power- 
law with no absorption, roughly consistent with the cosmic 
X-ray background, is appropriate. 

We can then use this new identification technique to 
make inferences about the heavily obscured AGN popula- 
tion from currently available X-ray survey data. Here we 
apply the scheme to the Chandra Deep Field North survey 
(CDFN). [Alexander et al.| ( |2003| ) present photometric data 
from this survey in the 0.5-1 keV, 1-2 keV and 2-8 keV bands. 
These correspond exactly to our band selection for z = 1 
sources. Redshift data for this catalogue have been provided 



by Trouille et al. (2008), including photometric redshifts. 
Using these, we select sources with 0.8 < z < 1.2. The hard- 
ness ratio plot for this survey is presented in Fig. [3] Here we 
also show tracks of a simple absorbed power-law, calculated 
using the zwabs*power-law models in xspec with r=1.7, 1.9 
and 2.1, from A r H =10 20 -10 24 cm" 2 . 

We find that for the 107 0.8 < z < 1.2 CDFN sources we 
have selected, 65 (61%) lie above the selection line. We can 
also use these hardness ratios to infer information about the 
spectral complexity present in these sources. Only 39 (36%) 
of the CDFN sources we study here show spectra consistent 
with a simple absorbed power-law, where the majority of 
sources present more complex spectra, due to soft excess 
emission. 



1 an online calculator can be found here 
https:/ /heasarc. gsfc.nasa.gov/Tools/w3pimms. html 



4 DISCUSSION AND SUMMARY 

In this paper we have aimed to identify an effective colour- 
colour selection for obscured AGN from X-ray data and im- 
prove upon the use of the classical hardness ratio. In Fig. 
[4] we plot histograms of the classical HR (with photon flux, 
rather than counts) for unabsorbed sources, absorbed Comp- 
ton thin sources and Compton thick sources from |Brightman| 



& Nandra ( 2011a I. This confirms our view that the hardness 



ratio is not an effective measure of obscuration, as all Comp- 
ton thick sources appear in fact soft, while both Compton 
thin and unabsorbed sources show the full range of hardness 
ratio. 

We have shown that colour-colour selection is able to 
break these degeneracies. The choice of bands for our colour- 
colour selection stems from the observation that while ob- 
scured AGN do appear hard above ~4 keV, they are usually 
dominated by a soft excess component below 4 keV. This 
component is often scattered nuclear light, which is believed 
to have undergone Thompson scattering off hot electrons 
far from the nucleus, and hence out-with the obscuring ma- 
terial. It typically represents only a few percent of the in- 
trinsic emission ( Turner et al.|1997 l. Soft X-ray emission in 
obscured AGN can also be attributed to circumnuclear gas 
photoionised by the AGN ( Guainazzi fc Bianchi|2007 1 , or to 
star-formation related processes in the host galaxy, such as 
X-ray binaries. In all cases, this emission is assumed to be 
not subjected to the obscuration of AGN, and thus appears 
soft in HR1. 

|Della Ceca et al.| ([l999) presented an analysis of 
hard X-ray selected sources from ASCA utilising very 
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Figure 1. An example set of six of our reference model spectra, 
for UGC 545, NGC 5506, NGC 262, NGC 2655, NGC 17 and 
NGC 3079 (top left to bottom right). These show a range in N-r 
values (logio(A r u/cm — 2 ) printed in the corner). The dashed lines 
show the intrinsic source component, and the dotted lines mark 
the boundaries of the bands that we use for the calculation of 
hardness ratios. 
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Figure 3. Our hardness ratio diagram applied to data from the 
CDFN. The solid red line is the selection criterion for Nn> 10 23 
cm -2 sources, described in the text. The solid lines black lines 
show the tracks of an absorbed power-law with (top to bottom) 



r=1.7, 1.9 and 2. 
left to right. 



1, from JVh=10 -10 cm , progressing from 




0.0 
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Figure 4. Histogram showing the distribution of the classical 
hardness ratio, using the soft band (0.5-2 keV) and the hard band 
(2-10 keV) photon fluxes from our sample of reference spectra. 
The blue histogram shows all sources, the green histogram shows 
A r jj> 10 23 cm -2 sources and the red histogram shows Compton 
thick sources. 



similar hardness ratios to ours. They also use these 
to infer absorption and spectral complexity in their 
sample. We have improved upon this technique to 
provide specific selection criteria for AGN with Nh> 
10 23 cm -2 as well as Compton thick AGN, and to 
extend the analysis to higher redshifts. Furthermore 



Noguchi et al. ( 2009 1 use similar hardness ratios in their 



study of the scattered fraction in AGN, also using data from 
2XMM. They use the hardness ratios to select AGN with 
very small scattered fractions, in the search for buried AGN 
first described in |Uedaet al.| ( |200"7| ). They also note, however, 
that heavily obscured AGN can be selected well using these 
hardness ratios, though this is not the focus of their work. 

Our results here provide a valuable technique for the 
identification of obscured AGN from photometry alone, 
when a spectral fit is not possible or badly constrained, and 
provides a direct method for investigating the nature of ob- 



scuration in AGN. As we have shown, it has applications for 
deep field, high redshift studies, where there are few good 
quality spectra, though above z~l, the 1-2 keV restframe 
band begins to be redshifted out of most X-ray telescope 
bandpasses. The redshift of the source should be known also, 
for correct application of the narrow bands. As Noguchi ct al. 
( |2009[ > suggested, it is well suited for application to wider 
and shallower surveys, such as the XMM-Newton serendipi- 
tous source catalogue, especially as the photometry provided 
in that catalogue matches well to our choice of bands. While 
this method does very well at separating sources 
with A?h> 10 23 cm -2 from sources with lower A^h, 
only X-ray spectroscopy can identify Compton thick 
sources, through the detection of the Fe Ka line. 
This will be made possible for the low fluxes consid- 
ered here with the high throughput X-ray spectral 
capabilities of ATHENA 
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We have also used our hardness ratio analysis to make 
inferences about the spectral complexity CDFN sources as 



Delia Ceca et al.| ( |l999| ) did for ASCA sources. We find 
that a significant number of CDFN sources at z ~1 show 
spectral shapes more complicated than a simple absorbed 
power-law due to soft-excess emission. This has implications 
for the X-ray selection of heavily obscured sources, as it 
shows that they can be detected in the soft band through 
their soft-excess emission. This is good news for the eRosita 
instrument will conduct an all-sky imaging survey up to 10 
keV. One of the primary science goals is to detect systemat- 
ically all obscured accreting Black Holes in nearby galaxies 
and many (up to 3 million) new, distant active galactic nu- 
clei. This technique provides a reliable tool for this purpose, 
especially if it is sensitive down to 0.2 keV, as this technique 
could be applied to sources out to z ~ 4. 
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